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Sequence Design to Minimize the Weighted 
Integrated and Peak Sidelobe Levels 
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Abstract —Sequences with low aperiodic autocorrelation side- 
lobes are well known to have extensive applications in active 
sensing and communication systems. In this paper, we consider 
the problem of minimizing the weighted integrated sidelobe 
level (WISL), which can be used to design sequences with 
impulse-like autocorrelation and zero (or low) correlation zone. 
Two algorithms based on the general majorization-minimization 
method are developed to tackle the WISL minimization problem 
and the convergence to a stationary point is guaranteed. In 
addition, the proposed algorithms can be implemented via fast 
Fourier transform (FFT) operations and thus are computationally 
efficient, and an acceleration scheme has been considered to 
further accelerate the algorithms. Moreover, the proposed meth¬ 
ods are extended to optimize the /',,-norm of the autocorrelation 
sidelobes, which lead to a way to minimize the peak sidelobe level 
(PSL) criterion. Numerical experiments show that the proposed 
algorithms can efficiently generate sequences with virtually zero 
autocorrelation sidelobes in a specified lag interval and can also 
produce very long sequences with much smaller PSL compared 
with some well known analytical sequences. 

Index Terms —Autocorrelation, majorization-minimization, 
peak sidelobe level, weighted integrated sidelobe level, unit- 
modulus sequences. 

I. Introduction 

S EQUENCES with good autocorrelation properties lie at 
the heart of many active sensing and communication 
systems. Important applications include synchronization of 
digital communication systems (e.g., GPS receivers or CDMA 
cellular systems), pilot sequences for channel estimation, 
coded sonar and radar systems and even cryptography for 
secure systems |TS-@|. In practice, due to the limitations 
of sequence generation hardware components (such as the 
maximum signal amplitude clip of analog-to-digital converters 
and power amplifiers), unit-modulus sequences (also known as 
polyphase sequences) are of special interest because of their 
maximum energy efficiency 13 - 

Let {x n }n =1 denote a complex unit-modulus sequence of 
length N, then the aperiodic autocorrelations of {a;„}^ =1 are 
defined as 

N-k 

Tk = X x *n x n+k = r*_ k , k = 0,..., N - 1. (1) 

n =1 

The problem of sequence design for good autocorrelation 
properties usually arises when small autocorrelation sidelobes 
(i.e., k ^ 0) are required. To measure the goodness of the 
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autocorrelation property of a sequence, two commonly used 
metrics are the integrated sidelobe level (ISL) 

JV-l 

ISL=^|r fe | 2 , (2) 

k= 1 

and the peak sidelobe level (PSL) 

PSL = maxdrfel}^ 1 . (3) 

Owing to the practical importance of sequences with low 
autocorrelation sidelobes, a lot of effort has been devoted to 
identifying such sequences. Binary Barker sequences, with 
their peak sidelobe level (PSL) no greater than 1, are perhaps 
the most well-known such sequences However, it is 
generally accepted that they do not exist for lengths greater 
than 13. In 1965, Golomb and Scholtz 0 started to investigate 
more general sequences called generalized Barker sequences, 
which obey the same PSL maximum, but may have complex 
(polyphase) elements. Since then, a lot of work has been done 
to extend the list of polyphase Barker sequences mui, and 
the longest one ever found is of length 77. It is still unknown 
whether there exist longer polyphase Barker sequences. Apart 
from searching for longer polyphase Barker sequences, some 
families of polyphase sequences with good autocorrelation 
properties that can be constructed in closed-form have also 
been proposed in the literature, such as the Frank sequences 
m, the Chu sequences lfj~4l . and the Golomb sequences fl5l . 
It has been shown that the PSL’s of these sequences grow 
almost linearly with the square root of the length N of the 

sequences ma, urn 

In recent years, several optimization based approaches have 
been proposed to tackle sequence design problems, see m- 
lf22l . Among them, |[22l and ED proposed to design unit- 
modulus sequences with low autocorrelation by directly mini¬ 
mizing the true ISL metric or a simpler criterion that is “almost 
equivalent” to the ISL metric. Efficient algorithms based on 
fast Fourier transform (FFT) operations were developed and 
shown to be capable of producing very long sequences (of 
length 10 4 or even larger) with much lower autocorrelation 
compared with the Frank sequences and Golomb sequences. 
Why the ISL metric was chosen as the objective in the 
optimization approaches? It is probably because the ISL metric 
is more tractable compared with the PSL metric from an 
optimization point of view. But as in the definition of Barker 
sequences, PSL seems to be the preferred metric in many 
cases. So it is of particular interest to also develop efficient 
optimization algorithms that can minimize the PSL metric. 
Additionally, in some applications, instead of sequences with 
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impulse-like autocorrelation, zero correlation zone (ZCZ) se¬ 
quences (with zero correlations over a smaller range) would 
suffice [1231 , l24l . In |fl8l , an algorithm named WeCAN 
(weighted cyclic algorithm new) was proposed to design 
sequences with zero or low correlation zone. 

In this paper, we consider the problem of minimizing the 
weighted ISL metric, which includes the ISL minimization 
problem as a special case and can be used to design zero or low 
correlation zone sequences by properly choosing the weights. 
Two efficient algorithms are developed based on the general 
majorization-minimization (MM) method by constructing two 
different majorization functions. The proposed algorithms can 
be implemented by means of FFT operations and are thus 
very efficient in practice. The convergence of the algorithms 
to a stationary point is proved. An acceleration scheme is 
also introduced to further accelerate the proposed algorithms. 
We also extend the proposed algorithms to minimize the t v - 
norm of the autocorrelation sidelobes. The resulting algorithm 
can be adopted to minimize the PSL metric of unit-modulus 
sequences. 

The remaining sections of the paper are organized as 
follows. In Section [III the problem formulation is presented. 
In Sections [HI] and II VI we first give a brief review of the MM 
method and then two MM algorithms are derived, followed by 
the convergence analysis and an acceleration scheme in Sec¬ 
tion [V] In Section [VI] the algorithms are extended to minimize 
the £p-norm of the autocorrelation sidelobes. Finally, Section 
I VIII presents some numerical results and the conclusions are 
given in Section IVIIII 

Notation: Boldface upper case letters denote matrices, bold¬ 
face lower case letters denote column vectors, and italics 
denote scalars. R and C denote the real field and the complex 
field, respectively. Re(-) and Im(-) denote the real and imag¬ 
inary part respectively, arg(-) denotes the phase of a complex 
number. The superscripts (-) T , (•)* and (•) H denote transpose, 
complex conjugate, and conjugate transpose, respectively, o 
denotes the Hadamard product. X, hJ denotes the (/-th, j- th) 
element of matrix X and Xi denotes the /-th element of vector 
x. X, ; denotes the /-th row of matrix X, X : J denotes the j- 
th column of matrix X, and denotes the submatrix of 

X from X l k to Xjj. Tr(-) denotes the trace of a matrix. 
diag(X) is a column vector consisting of all the diagonal 
elements of X. Diag(x) is a diagonal matrix formed with x 
as its principal diagonal. vec(X) is a column vector consisting 
of all the columns of X stacked. I„ denotes annxn identity 
matrix. 

II. Problem Formulation 

Let {iC-ri} rT= i denote the complex unit-modulus sequence to 
be designed and { f'k } k _ t ' be the aperiodic autocorrelations 
of {x n }n = i as defined in ([]}. then we define the weighted 
integrated sidelobe level (WISL) as 

N -1 

WISL= ^TtUfekl 2 , (4) 

fc =1 

where w k > 0, k = 1,..., N — 1. It is easy to see that the 
WISL metric includes the ISL metric as a special case by 
simply taking Wk = 1, k = 1,... ,N — 1. 


The problem of interest in this paper is the following WISL 
minimization problem: 

minimize WISL 

*■ (5) 

subject to | x n | = l,n = l,..., N, 

which includes the ISL minimization problem considered in 
ea as a special case and can be used to design zero (or low) 
correlation zone sequences by assigning larger weights for the 
sidelobes that we want to minimize. 

An algorithm named WeCAN was proposed in fT8l to tackle 
this problem. However, instead of directly minimizing the 
WISL metric as in ©, WeCAN tries to minimize an “al¬ 
most equivalent” criterion. Moreover, the WeCAN algorithm 
requires computing the square-root of an N x N matrix at the 
beginning and N FFT’s at each iteration, which could be costly 
for large N. In the next section, we will develop algorithms 
that directly minimize the original WISL metric, and at the 
same time are computationally much more efficient than the 
WeCAN algorithm. 

III. Sequence Design via 
Majorization-Minimization 

In this section, we first introduce the general majorization- 
minimization (MM) method briefly and then apply it to derive 
simple algorithms to solve the problem <[5}i. 

A. The MM Method 

The MM method refers to the majorization-minimization 
method, which is an approach to solve optimization problems 
that are too difficult to solve directly. The principle behind the 
MM method is to transform a difficult problem into a series 
of simple problems. Interested readers may refer to |j25l , |26l 
and references therein for more details (recent generalizations 
include |ZQ, l28l l. 

Suppose we want to minimize /(x) over X C C". 
Instead of minimizing the cost function /(x) directly, the 
MM approach optimizes a sequence of approximate objective 
functions that majorize /(x). More specifically, starting from a 
feasible point x^°\ the algorithm produces a sequence {x^ fc ^} 
according to the following update rule: 

x (fc+i) g argmin u(x, x^), (6) 

xeA 

where x !k> is the point generated by the algorithm at iteration 
k, and u(x, x^ fc ^) is the majorization function of /(x) at 
x ( fc ). Formally, the function u(x, x^ fc ^) is said to majorize the 
function /(x) at the point if 

w(x,x (fe) ) > /(x), Vx G X, (7) 

u(xW,x(‘)) = /(x (fe )). (8) 

In other words, function M(x,X< fe )) is an upper bound of /(x) 
over X and coincides with /(x) at x tk >. 

To summarize, to minimize /(x) over X C C ra , the main 
steps of the majorization-minimization scheme are 

1) Find a feasible point x® and set k = 0. 

2) Construct a function it(x, x^) that majorizes /(x) at 

#). 
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3) Let x( fc+1 ) g argmin u(x, x( fc )). 

xe* 

4) If some convergence criterion is met, exit; otherwise, set 
k = k + 1 and go to step (2). 

It is easy to show that with this scheme, the objective value 
is monotonically decreasing (nonincreasing) at every iteration, 
i.e., 

J( x ( fc +!)) < u(x^ k+1 \x^) < M(xW,x (i) ) = /(x (fc) ). (9) 

The first inequality and the third equality follow from the the 
properties of the majorization function, namely © and © 
respectively and the second inequality follows from ©. The 
monotonicity makes MM algorithms very stable in practice. 

B. WISL Minimization via MM 


and denote the objective function of (ITTI i by /(X), then 
/(X) = vec(X) i/ Lvec(X), which is clearly a quadratic 
function in X. According to Lemma [U we can construct a 
majorization function of /(X) by simply choosing a matrix 
M such that M ^ L. A simple choice of M can be 
M = A max (L)I, where A max (L) is the maximum eigenvalue 
of L. Owing to the special structure of L, it can be shown 
that A max (L) can be computed efficiently in closed form. 

Lemma 2. Let L be the matrix defined in <m. Then the max¬ 
imum eigenvalue ofL is given by A max (L) = ma Xk{w k {N — 
k)\k = l,...,N-l}. 

Proof: It is easy to see that the set of vectors 
are mutually orthogonal. For k = 1 — 
N,... ,N — 1, we have 


To solve the problem © via majorization-minimization, the 
key step is to find a majorization function of the objective such 
that the majorized problem is easy to solve. For that purpose 
we first present a simple result that will be useful later. 

Lemma 1 ||22| . Let L be an n x n Hermitian matrix and M 
be another nxn Hermitian matrix such that M >z L. Then for 
any point xo G C n , the quadratic function x fl Lx is majorized 
by x fl Mx + 2Re (x H (L — M)xo) + x^(M — L)xo at xo- 

Let us define Ufc, k = 0, — 1 to be N x N 

Toeplitz matrices with the fcth diagonal elements being 1 and 
0 elsewhere.Noticing that 

r k = Tr(U fc xx H ), k = 0,... ,N - 1, (10) 

we can rewrite the problem © as 


N -1 


JV-1 

minimize > 

X.x ' 

k =l 


w fe |TA(U fe X)| 2 


subject to X = xx fl 

\x n \ = l,n=l,...,N. 


( 11 ) 


Let U_fc = ~Uj. k = 1,...,, N — 1, it is easy to see that 
Tr(U_/;X) = r_fc = r* k , so we can rewrite the problem ([TH 
in a more symmetric way: 


minimize 

X.x 


Wfc |Tr(U fc X)| 2 


N -1 

I E 

k=l — N 

subject to X = xx H 

|x„| = 1, n = 1,..., TV, 


( 12 ) 


Lvec(Ufc) = ^ tu j ve c (U j) ve c (U j) ^ vec (U ^) 

j=1 ~ N (15) 

= uifcvec(U fe )vec(U fe ) ff vec(U fe ) 

= w k (N - |fc|)vec(U fe ). 

Thus Wk(N — |fc|), k = 1 — TV,..., N — 1 are the nonzero 
eigenvalues of L with corresponding eigenvectors vec(Ufc), 
k = 1 — N,..., N — 1. Since W-k = Wk and wq = 0, the 
maximum eigenvalue of L is given by maxfc{rUfe(7V — k)\k = 

1,...,JV-1}. ■ 

Then given X^- 1 = x^(x^) H at iteration l, by choosing 
M = A max (L)I in Lemma Q] we know that the objective of 
(fl3l) is majorized by the following function at X !, A 

U1 (X,XW) 

= A max (L)vec(X)"vec(X) 

+ 2Re(vec(X) ff (L - A max (L)I)vec(X«)) 

+ vec(X(')) H (A max (L)I - L)vec(xW). 

Since vec(X) ff vec(X) = (x H x) 2 = TV 2 , the first term of 
(fT6l) is just a constant. After ignoring the constant terms, the 
majorized problem of | |T3| > is given by 

minimize Re(vec(X) ff (L — A max (L)I)vec(X^)) 

subject to X = xx H 

\x n \ = 1, n= 1,...,TV. 

(17) 

Substituting L in ( 1 1 4| i back into (fTTI i. the problems becomes 


(16) 


N-l 


where w-k = Wk and wq = 0. Since Tr(UfcX) = 
vec(X) ff vec(Ufc), the problem (fl2l i can be further rewritten 
using vec notation as (the constant factor A is ignored): 


m 


inimize ^ w k Re(TV(U_ fc X('))Tr(U fc X)) 


-A 1 m I(L)Tr(X( ; )X) 

subject to X = xx H 


(18) 



N-l 


I 

X n \ = 

1, n= 1,...,TV, 

minimize 

iu fe vec(X) ff vec(Ufc)vec(Ufc)' f/ vecl 

;x) 



n\ _ 

x,X 

k—l — N 


Since Tr(U_ fc X^l) = 

r_ k , the problem (118l> 

subject to 

X = xx H 


as 




\x n \ = 1, n = 1,.. .,TV. 




( ( N - 1 



(13) 

minimize 

Re | 

Tr E iv k r«l U fc 

Let us define 


x,X 


\ \k=l-N 


N-l 



-A max (L)Tr(X«X) 


L = ^2 ^/evec(Ufc)vec(U/ c ) if 

(14) 

subject to 

X = 

= xx w 


k=l—N 



| %n | 

= 1, n = 


( 19 ) 
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which can be further simplified as 


minimize x H (R — A max (L)x® (x^) H ) x 

X V 7 

subject to \x n \ = 1 , n = 1 , ..., 

where 

N -1 

R= ]T Wfcr^Ufc 

fc=l-iV 

0 

tuir^ 0 

(0 (0 

_'Wn-i1'n-i ■ ■ ■ w i r i 

is a Hermitian Toeplitz matrix, and are t ^ le 

autocorrelations of the sequence {x^/ }^ =1 . 

It is clear that the objective function in ( l20t is quadratic 
in x, but the problem (l20l > is still hard to solve directly. So 
we propose to majorize the objective function of problem (|2(H > 
at x® again to further simplify the problem that we need to 
solve at each iteration. Similarly, to construct a majorization 
function of the objective, we need to find a matrix M such 
that M y R — A ma x(L)xW(xW) ff . As one choice, one 
may choose M = A max (R - A max (L) x (0 (x.^) H ) I, as in 
the first majorization step. But in this case, to compute the 
maximum eigenvalue of the matrix R — A max (L)x(')(x( i )) i/ , 
some iterative algorithms are needed, in contrast to the simple 
closed form expression in the first majorization step. To 
maintain the simplicity and the computational efficiency of 
the algorithm, here we propose to use some upper bound of 
Amax (R ^ A ma x(L)x® (x^) ff ) that can be easily computed 
instead. To derive the upper bound, we first introduce a useful 
result regarding the bounds of the extreme eigenvalues of 
Hermitian Toeplitz matrices fl29l . 


(0 ' 
WN-ir{_ N 


(0 

W\r_\ 


( 21 ) 




Lemma 3. Let T be an N x N Hermitian Toeplitz matrix 
defined by {tk}^=o as follows 


'to t\ ... t* N _, 

ti t 0 ■ 

: ■ . ■ . 

1 i 

fjv-1 • ■ • t\ to 


and F be a 2 N x 2TV FFT matrix with F m , n = e , 0 < 
m, n < 2N. Let c = [fo, ti, • • • , ijv-i, 0, t’^_ 1 , • • • , t\] T and 
fi = Fc be the discrete Fourier transform of c. Then 


Amax(T) < 


1 

2 


max u 2i + max u 2 i-i 
\l<i<N l<i<N 


( 22 ) 


Amin(T) > 


1 

2 


( min p 2 i + min p 2 i -1 

\l<i<N l<i<N 


(23) 


Proof: See (29]. ■ 

Since the matrix R is Hermitian Toeplitz, according to 
Lemma [3j we know that 

A max (R) < x ( max p 2i + max /z 2 i-i) , (24) 

2 Vi<j<iv 1 <i<N J 


where //, = Fc and 

c = [0, wir[ l \ ..., 0, lojv-irfljv) ■ • ■, wir ( '^ 1 ] T . 

(25) 

Let us denote the upper bound of A max (R) at the right hand 
side of (l24l > by A u , i.e., 

A,, = - f max ji 2i + max p 2 i-i) . (26) 

2 \l<i<N l<i<N J 


Since A max (L) > 0, it is easy to see that 

A„ > Amax (R) > Amax (R - A nl ax (L)x« (x« ) H ) . (27) 

Thus, we may choose M = A„I in Lemma|]]and the objective 
of (I20l > is majorized by 

« 2 (x,xW) 

= A u x H x+2Re ^x' ff (R-A m ax(L)x (n (x (,) } ff -A„I)x (i) ^ 
+ (x«) ff (A U I - R+A max (L)x«(xW) ff )x«. 

(28) 

Since x H x = N, the first term of (f28l> is a constant. Again 
by ignoring the constant terms, we have the majorized problem 
of ®: 

minimize Re (x H (R — A max (L)x(^(x® — A„I)x^)) 

X V 7 

subject to \x n \ = 1, n = 1,..., N, 

(29) 

which can be rewritten as 


minimize llx — y|L 

x 1 J (3Q) 

subject to \x n \ = 1, n = 1,..., N, 

where 

y = -(R - A max (L)x (n (x ( ' ) ) i/ - A u I)x (i) 

= (A max (L) 7 V + A u )x« ^Rx«. ( 31 ) 

It is easy to see that the problem (l30t has a closed form 
solution, which is given by 

x n =e jais ^ yn \ n=l,...,N. ( 32 ) 


Note that although we have applied the majorization- 
minimization scheme twice at the point x^', it can be viewed 
as directly majorizing the objective function of (0 at x® by 
the following function: 

«(x,xW) 

N—l 

= w 2 (x,x (z) ) + A max (L)iV 2 - Yw k \rP\* 

k=1 

N—l 

= - 2Re (x*y) +2iV(A max (L)7V + A u )-3^^|r«| 2 , 

k =1 

(33) 

and the minimizer of u(x. x^)) over the constraint set is given 
by (EJ. 

According to the steps of the majorization-minimization 
scheme described in section IIII-AI we can now readily have 
a straightforward implementation of the algorithm, which at 
each iteration computes y according to (PH and update x via 
(132L It is easy to see that the main cost is the computation of y 
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in (ED- T° obtain an efficient implementation of the algorithm, 
here we further explore the special structures of the matrices 
involved in the computation of y. 

We first notice that to compute X Ul we need to compute 
the FFT of the vector c in ( l25l > and the autocorrelations 
{ r k ) }k=i-N °f {^n^ \n=i are needed to form the vector c. 
It is well known that the autocorrelations can be computed 
efficiently via FFT (IFFT) operations, i.e., 


r r (0 r (0 
L'o )'i > 




„(0 


1 

2 N 


F h F[x® t , 0 1x jv 


W-l> r l-jV> • ■ • > ' -1 
T 2 


(OlT 


(34) 


where F is the 2N x 2N FFT matrix and |-| 2 denotes the 
element-wise absolute-squared value. Next we present another 
simple result regarding Hermitian Toeplitz matrices that can 
be used to compute the matrix vector multiplication Rx ,!) 
efficiently via FFT (IFFT). 


Lemma 4. Let T be an N x N Hermitian Toeplitz matrix 
defined as follows 

't 0 t\ ... t* N _f 

rj, ti to '• : 

: ■. •. 4.* 

tN-l ■ ■ ■ t\ to 

and F be a 2N x 2N FFT matrix with F mjn = 
e~ 3 2 n ; 0 < m,n < 2N. Then T can be decom¬ 
posed as T = 2 ^F ff 1 . Ar Diag(Fc)F : i : jv, where c = 

[■to,tiAV- 


Proof: The N x N Hermitian Toeplitz matrix T can be 
embedded in a circulant matrix C of dimension 2 N x 2 N as 
follows: 


C = 


T W 
W T ’ 


(35) 


where 


0 tjv-i 


t 


l 


W = 


t 


* 

N— 1 


t 


* 

1 


t 


* 

N—l 


tN-l 

0 


(36) 


The circulant matrix C can be diagonalized by the FFT matrix 

ED, i-e., 

C = :4F if Diag(Fc)F 1 (37) 


where c is the first column of C, i.e., c = 
[to,ti, ■ ■ ■ , t^-i, Since the matrix T 

is just the upper left N x N block of C, we can easily obtain 

T = 277 F ^l:iV Dia g( Fc ) F :,l:Af. ■ 

Since the matrix R is Hermitian Toeplitz, from Lemma [4] 
we easily have 

R = 2F F fl:iV Di ag( Fc ) F :,l:iV, (38) 

where c is the same as the one defined in (l25l >. which can be 
reused here. With the decomposition of R given in ( 1381 ). it is 
easy to see that the matrix vector multiplication Rx® can be 
performed by means of FFT (IFFT) operations. 


Now we are ready to summarize the overall algorithm 
and it is given in Algorithm Q] The algorithm falls into the 
general framework of MM algorithms and thus preserves the 
monotonicity of such algorithms. It is also worth noting that 
the per iteration computation of the algorithm is dominated 
by four FFT (IFFT) operations and thus computationally very 
efficient. 


Algorithm 1 MWISL - Monotonic minimizer for Weighted 
ISL. 


Require: sequence length N, weights {tvt > 
l: Set l = 0, initialize x(°). 

2: Xl = maxfc{wfc(N - k)\k = 1,... ,N - 1} 

3: repeat 

4: f = F[x« T ,0i xW ] T 

5: r = 2jF ff |f| 2 

6: C = ro [0, WI, . . .,W N -1,0,W N -1, . . .,Wi] T 

7: pL = Fc 

8: A„ = max u 2 j + max u 2 i-i) 

z X 1 <i<N l<i<N ' 


(0 - 


y = x 

„( ; +l) _ j arg 
— C- 




2N(\ l N+\ u ) 

l) , n = 1,..., N 

l 4-l + l 

until convergence 


IV. WISL Minimization with an Improved 
Majorization Function 

As described in the previous section, the proposed algorithm 
is based on the majorization-minimization principle, and the 
nature of the majorization functions usually dictate the per¬ 
formance of the algorithm. In this section, we further explore 
the special structure of the problem and construct a different 
majorization function. 

Notice that to obtain the simple algorithm in the previous 
section, a key point is that the first term of the majorization 
function ui(X,X(‘>) in © is just a constant and can be 
ignored, which removes the higher order term in the objective 
of the majorized problem (ITTb . In the previous section, we have 
chosen M = A max (L)I + L such that vec(X) ff Mvec(X) = 
Amax(L)vec(X) ff vec(X) is a constant over the constraint set. 
But it is easy to see that, to ensure the term vec(X) ff Mvec(X) 
being constant, it is sufficient to require M being diagonal, i.e., 
choosing M = Diag(b) >: L. To construct a tight majorization 
function, here we consider the following problem 

minimize Tr(Diae(b) — L) 

b (39) 

subject to Diag(b) + L, 

i.e., we choose the diagonal matrix Diag(b) that minimizes 
the sum of eigenvalues of the difference Diag(b) — L. Since 
L is a constant matrix, the problem ( l39l > can be written as 

minimize b T l„ 

b (40) 

subject to Diag(b) + L. 

The problem ( l4Qt is an SDP (semidefinite programming), and 
there is no closed form solution in general. However, due to 
the special properties (i.e., symmetry and nonnegativity) of the 
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matrix L in (IT4l) . it can be shown that the problem ( l40l ) can be 
solved in closed form and it is given in the following lemma. 

Lemma 5. Let L be an n x n real symmetric nonnegative 
matrix. Then the problem 

minimize b 1 l n 

b (41) 

subject to Diag(b) L 

admits the following optimal solution: 

b*=Ll n . (42) 

Proof: Since Diag(b) - L 0, we have 

x T (Diag(b) — L) x > 0, Vx. 


where mat(-) is the inverse operation of vec(-), following the 
derivation in the previous section, we can simplify ( |45| > to 

minimize x fl (R-Bo(x^(x®) ff ))x 

X v \ \ j // (47 ) 

subject to \x n \ = 1, n = 1,. .., AT, 
where R is defined in (EH and 

B = mat(b) 

= mat(Ll) 

/ JV-i \ 

= mat I 55 w’fcvec(Ufc)vec(U fc ) H l J 

\k=l-N ) 

N -1 

= MN-\k\)V k . 


By choosing x = 1„, we get b 2 1„ > l^Ll n , with the 
equality achieved by b = Ll n . So it remains to show that 
b = Ll n is feasible, i.e., Diag(Ll n ) >z L. Since for any 
x G R n , we have 

x T (Diag(Ll„) — L) x 

= y ( x i 55 Lij)— 55 xiLijXj 

i j i,j 

= 55 — XiLijXj) 


hi 





> 0 , 


where the third equality follows from the symmetry of L and 
the last inequality follows from the fact that L is nonnegative, 
the proof is complete. ■ 

Then given = xOl(xO)) 22 at iteration l, by choosing 
M = Diag(b) = Diag(Ll) in LemmaQ] the objective of (ITTl i 
is majorized by the following function at X 1- ^: 

ui(X,X w ) 

= vec(X) 22 Diag(b)vec(X) 

+ 2Re(vec(X) ff (L — Diag(b))vec(X^^)) 

+ vec(xW) ff (Diag(b) - L)vec(X®). 

As discussed above, the first term of (l44l > is a constant and by 
ignoring the constant terms we have the majorized problem 
given as follows: 

minimize Re(vec(X) ff (L — Diag(b))vec(X®)) 

subject to X = xx 22 (45) 

\x n \ = l,n=l,...,N. 

Since 

Re(vec(X) 22 Diag(b)vec(X^- ) )) 

= Re(vec(X) 22 (b o vec(X^^))) 

= Tr(xx 22 mat(b o vec(x^(x^) 22 ))) ^6) 

= x 22 ^mat(b) o (xW(xW) 22 )) x, 


The objective in (l47l) is quadratic in x and similar as 
before we would like to find a matrix M such that 
M A R - B o (xW(xW) 22 ). For the same reason as in 
the previous section, here we use some upper bound of 
Amax (R — B o (x^^x®) 22 )) to construct the matrix M. 
We first present two results that will be used to derive the 
upper bound. The first result indicates some relations between 
Hadamard and conventional matrix multiplication ED- 

Lemma 6. Let A, Be Qmxn an( j x g p/ len p le 
diagonal entry of the matrix ADiag(x)B 2 coincides with the 
ith entry of the vector (A o B)x, i = 1,..., m, i.e., 

(A o B)x = diag (ADiag(x)B T ) . (48) 

Then the second result follows, which reveals a fact regard¬ 
ing the eigenvalues of the matrix B o (x^^x^)) 22 ). 

Lemma 7. Let B be an N x N matrix and x €E C Y with 
\x n \ = 1, n = 1,..., N. Then B o (xx 22 ) and B share the 
same set of eigenvalues. 


Proof: Suppose A is an eigenvalue of B and z is the 
corresponding eigenvector, i.e., Bz = Az, then 

(B o (xx 22 )) (x o z) 

= diag(BDiag(x o z)(xx 22 ) T ) 

= diag(B(xozox*)x T ) 

= diag(Bzx T ) 

= diag(Azx T ) 

= A(x o z), 

which means A is also an eigenvalue of the matrix Bo (xx 22 ), 
with the corresponding eigenvector given by (x o z). ■ 

With Lemma [7] we have 


A 

^max 
= Amax 


^Bo(x (l) (x (,) ) 22 
(R) - A min (B O ( x ^^ x '^) 22 
(R) - Amin (B) . 


(50) 


Noticing that the matrix B is symmetric Toeplitz, according 
to Lemma [3] we know that 


An 


a( B ) > \ 


min u 2 i + min v 2 i-i 

l<i<N l<i<N 


(51) 
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where v = Fw and 


w = [0, w\[N - 1 ),... ,w N - 1 ,0,w N - 1 ,.. .,W!(N - 1)] T . 

(52) 

By defining 

As = 7: ( min is 2i + min v 2 i- 1 ) , (53) 

2 \l<i<lV 1 <i<N J 

we now have 

A max (R-Bo( x W(x«) ff )) <\ u -\b, (54) 

where A„ is defined in (l26t and by choosing M = (A u — As)I 
in Lemma [I] we know that the objective of ( 147 b is majorized 
by 


« 2 (x,x«) 

= (A„ - A B )x ff x 

+ 2Re (x ff (R B o (x^(x^) H ) -(A u -A b )I)xW) 
+ (x«) ff ((A u - A b )I - R+B o ( x (0(xW)H)) x (0. 


Algorithm 2 MWISL-Diag - Monotonic 
Weighted ISL. 

minimizer for 

Require: sequence length N, weights {wt > 0} fc _j 1 

1 

Set l = 0, initialize x^ 0 -*. 


2 

w= [0, wi(N-l),... ,wn-i, 0, Wjv-i, 


3 

v — Fw 


4 

P = F 5: n{ VO ( F :,l:Vl)) 


5 

A B = b (mini<j<jv u 2i + mini<j<jv V2i- 

-i) 

6 

repeat 


7 

f = F[x« t ,0 1xJV ] t 


8 

r=^F H \f\ 2 


9 

c = ro[0, wi,..., wn-i,0, iojv- 1 , ■ • ■ 

,wi] T 

10 

fi = Fc 


11 

A„ = i( max jx 2i + max /r 2 j-i) 



2 y l<i<N l<i<N ' 


12 

„m , P ox(0 '- F "i : Ar(r‘°f) 
y _ x ^ 2N(X u —X b ) 


13 

x { n +1) n = l,...,JV 


14 

l <-l + 1 


15 

until convergence 



By ignoring the constant terms, the majorized problem of (l47l > V. Convergence Analysis and Acceleration 
is given by SCHEME 


minimize Re [x 11 (R—B o (x^(x^) ff ) — (A u — As)I)x(^) 
subject to \x n \ = 1, n = 1,..., N. 

(55) 

Similar to the problem ( l29l > in the previous section, ( l55l > also 
admits a closed form solution given by 


n = 1,-TV, 


(56) 


where 


y = (A u —A s )x (i) +(Bo (x (1) (x (,) ) ff ))x (1) -Rx (l >. 

(57) 

It is worth noting that the y in ( 1571 ) can be computed efficiently 
via FFT operations. To see that, we first note that Rx-7 can 
be computed by means of FFT as described in the previous 
section. According to Lemma [6] we have 

(Bo( x (‘)(x«) fl ))x (i) 

= diag(BDiag(x«) (x® (x^) H ) T ) 

= diag(B(x«o (x«)*)(x«) T ) (58) 

= diag(Bl(x^) T ) 

= (Bl)ox®. 


A. Convergence Analysis 

The MWISL and MWISL-Diag algorithms given in Al¬ 
gorithm Q] and [2] are based on the general majorization- 
minimization framework, thus according to subsection IIII-AI 
we know that the sequence of objective values (i.e., weighted 
ISL) evaluated at {x^} generated by the algorithms is non¬ 
increasing. And it is easy to see that the weighted ISL metric 
in (Q} is bounded below by 0, thus the sequence of objective 
values is guaranteed to converge to a finite value. 

Now we further analyze the convergence property of the 
sequence {x^)} itself. In the following, we will focus on the 
sequence generated by the MWISL algorithm (i.e. Algorithm 
0}. and prove the convergence to a stationary point. The same 
result can be proved for the MWISL-Diag algorithm (i.e. 
Algorithm [3 similarly. 

To make it clear what is a stationary point in our case, we 
first introduce a first-order optimality condition for minimizing 
a smooth function over an arbitrary constraint set, which 
follows from |[32l . 

Proposition 8. Let f : R" -A If be a smooth function, and 
let x* be a local minimum of f over a subset X of R". Then 

V/(x*) T z > 0, Vz e Tx(x*), (60) 


Since B is symmetric Toeplitz, by Lemma |4] we can decom¬ 
pose B as 


where Txix*) denotes the tangent cone of X at x*. 


B = ^ F ^i : Ar Di ag(FW)F :i i :JV , (59) 

and thus B1 can also be computed efficiently via FFT opera¬ 
tions. We further note that we only need to compute B1 once. 
The overall algorithm is then summarized in Algorithm |2] for 
which the main computation of each iteration is still just four 
FFT (IFFT) operations and thus of order 0(N log N). 


A point x £ X is said to be a stationary point of the problem 

minimize fix) 
xex J v 

\i it satisfies the first-order optimality condition (f60t . 

To facilitate the analysis, we further note that upon defining 

x = [Re(x) T , Im(x) T ] T 


( 61 ) 










U *=2 


U fc 


u£ 


u fc 




u, = 


1 


0 

U^-U fc 


Ufc - u; 

0 


(62) 


(63) 


and based on the expression of /> in (ITOl) . it is straightforward 
to show that the complex WISL minimization problem ([5} is 
equivalent to the following real one: 


minimize 


JV-l 

E 

k= 

x 2 


]w k ^x T U fc x) + ^x T UfcX^ ^ 
= 1, n = 1,..., N. 


(64) 

K=1 

subject to x'i + x 2 n+N 
We are now ready to state the convergence properties of 
MWISL. 


Theorem 9. Let {x^} be the sequence generated by the 
MWISL algorithm in Algorithm [7] Then every limit point of 
the sequence {x^} is a stationary point of the problem <0. 

Proof: Denote the objective functions of the problem 0 
and its real equivalent (l64t by /(x) and /(x), respectively. 
Denote the constraint sets of the problem 0 and (l64l) by C and 
C, respectively, i.e., C = {x G C^l \x n \ = 1, n = 1,..., N} 
and C = {x G R 2Ar | x 2 n + x 2 n+N = 1, n = 1,..., N}. From 
the derivation of MWISL in subsection Illl-Bl we know that, 
at iteration l, /(x) is majorized by the function m(x,x(^) in 
(l33l > at x.' 1 ' 1 over C. Then according to the general MM scheme 
described in subsection IIII-AI we have 

/(x^ +1 )) < u(x l ' l+1 \x 1 ' 1 ' 1 ) < u(xS l \yS 1 ^) = /(x*^), 

which means {/(x^))} is a nonincreasing sequence. 

Since the sequence (x^)} is bounded, we know that it has 
at least one limit point. Consider a limit point x(°°) and a 
subsequence {x^A} that converges to x(°°\ we have 

u(x ( ^+ l) ,x ( ^+ l) ) = /(x ( ^+ l} ) < /(x^ +1 )) 

< < w(x, x^),Vx G C. 

Letting j —> +oo, we obtain 

ii(x (oo) ,x H )<«(x,x w ),VxeC, (65) 

i.e., x(°°) is a global minimizer of u(x, x(°°)) over C. With 
the definitions of x, Ufc and Ufc given in ( 16 1 1 . (1621 ) and < l63b . 
and by ignoring the constant terms in w(x, x^°°)), it is easy to 
show that minimizing it(x, x(°°)) over C is equivalent to the 
following real problem: 

minimize 4x T d — 2(A max (L)iV + A u )x T x( 0 °) 

* _ (66) 
subject to xgC, 

where x(°°) = [Re(x(°°() T ,Im(x(°°^) T ] T and 
JV-l 

d = ^ w k (x (oo)T UfcX (oo) Ufc + x (oo)T UfcX (oo) Ufc^ x (oo) . 

fc=l 

(67) 

Since x^ 00 ^ minimizes it(x, x^°°^) over C, x^°°^ is a global 
minimizer of ( l66t and since x T x is just a constant over C, 
x(°°) is also a global minimizer of 

minimize 4x T d — (A max (L)A^ + A„) (2x T x(°°( — x T x) 

subject to xgC. 

(68) 


Then as a necessary condition, we have 

V{t(ic^) T z > 0, Vz e T c ~(x(°°>), (69) 

where m(x) denotes the objective function of (l68l >. It is easy 
to check that 

V/(x (oo) ) = V5(x (oo >) = 4d. 

Thus we have 

V/(x(°°)) T z > 0, Vz G Tc(x (oo) ), (70) 

implying that x(°°) is a stationary point of the problem (l64l >. 
Due to the equivalence of problem (|64]> and 0, the proof is 
complete. ■ 

B. Acceleration Scheme 

In MM algorithms, the convergence speed is usually dictated 
by the nature of the majorization functions. Due to the 
successive majorization steps that we have carried out in the 
previous sections to construct the majorization functions, the 
convergence of MWISL and MWISL-Diag seems to be slow. 
In this subsection, we briefly introduce an acceleration scheme 
that can be applied to accelerate the proposed MM algorithms. 
It is the so called squared iterative method (SQUAREM), 
which was originally proposed in lf33ll to accelerate any Ex¬ 
pectation-Maximization (EM) algorithms. SQUAREM adapts 
the idea of the Cauchy-Barzilai-Borwein (CBB) method lf34ll . 
which combines the classical steepest descent method and the 
two-point step size gradient method [35], to solve the nonlinear 
fixed-point problem of EM. It only requires the EM updating 
scheme and can be readily implemented as an off-the-shelf 
accelerator. Since MM is a generalization of EM and the 
update rule is also just a fixed-point iteration, SQUAREM 
can be easily applied to MM algorithms after some minor 
modifications. 

Suppose we have derived an MM algorithm to minimize 
/(x) over X C C n and let Fmm(') denote the nonlinear 
fixed-point iteration map of the MM algorithm: 

x( fc+1 > =F MM (x (fc) ). (71) 

For example, the iteration map of the MWISL algorithm 
is given by (l32l >. Then the steps of the accelerated MM 
algorithm based on SQUAREM are given in Algorithm [3] 
A problem of the general SQUAREM is that it may vio¬ 
late the nonlinear constraints, so in Algorithm [3 we need 
to project wayward points back to the feasible region by 
Vx{-)- For the unit-modulus constraints in the problem under 
consideration, the projection can be done by simply applying 
the function ef aTg< ^ element-wise to the solution vectors. 
A second problem of SQUAREM is that it can violate the 
descent property of the original MM algorithm. To ensure 
the descent property, a strategy based on backtracking has 
been adopted in Algorithm [3] which repeatedly halves the 
distance between a and —l:a ■<— (a — l)/2 until the descent 
property is maintained. To see why this works, we first note 
that Vx (x( fc ) — 2ar + a 2 v) = X 2 if a = — 1. In addition, 
since /(x 2 ) < /(x^j due to the descent property of original 
MM steps, /(x) < /(x* fe )j is guaranteed to hold as a —> —1. 
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It is worth mentioning that, in practice, usually only a few 
back-tracking steps are needed to maintain the monotonicity 
of the algorithm. 


Algorithm 3 The acceleration scheme for MM algorithms. 
Require: parameters 
1: Set k = 0, initialize x*- 0 ^. 

2 : repeat 

3 : Xi = FmM (x (fc) ) 

4: X 2 = Fmm (xi) 

5: r = Xi — x( k ) 

6: v = x 2 — xi — r 

7: Compute the step-length a = 

8: x = Vx (x( fe ) — 2ar + ct 2 v) 

9: while /(x) > /(x^) do 

10: a <- (a - l)/2 

11: x = Vx (x( fc ) — 2ar + a 2 v) 

12 : end while 

13: x (fc+1) = X 

14: k £- k + 1 

15: until convergence 


VI. Minimizing the £ p -norm of Autocorrelation 
Sidelobes 


In previous sections, we have developed algorithms to 
minimize the weighted ISL metric of a unit-modulus sequence. 
It is clear that the (unweighted) ISL metric is just the squared 
1 2 -norm of the autocorrelation sidelobes and in this section 
we would like to consider the more general £ p -norm metric of 
the autocorrelation sidelobes defined as 

/N-i X 1 /? 5 

5 >r (72) 


with 2 < p < oo. The motivation is that by choosing different 
p values, we may get different metrics of particular interest. 
For instance, by choosing p -A +oo, the ( p -norm metric tends 
to the f'oo-norm of the autocorrelation sidelobes, which is 
known as the peak sidelobe level (PSL). So it is well motivated 
to consider the more general f p -norm (2 < p < oo) metric 
minimization problem 

( N-l 

Ei rfc i p 

fc =1 

subject to \x n \ = 1, n = 1,..., N, 
which is equivalent to 



JV-l 

minimize 

k =1 

subject to \x n \ = 1, n = 1,..., N. 


(74) 


It is easy to see that by choosing p = 2, problem (1741) becomes 
the ISL minimization problem. It is also worth noting that we 
can easily incorporate weights in the objective (i.e., weighted 
^ p -norm) as in (0]). 


To tackle the problem (f74l > via majorization-minimization, 
we need to construct a majorization function of the objective 
and the idea is to majorize each \r k \ p , k = 1,..., N — 1 by 
a quadratic function of |r*fe |. It is clear that when p > 2, 
it is impossible to construct a global quadratic majorization 
function of \r k \ p . However, we can still majorize it by a 
quadratic function locally based on the following lemma. 

Lemma 10. Let f(x) = x p with p > 2 and x £ [0, t]. Then 
for any given Xq £ [0, t), f(x) is majorized at Xq over the 
interval [0, t] by the following quadratic function 

ax 2 + (pXq 1 — 2axo)x + ax g — (p — 1):£q, (75) 


where 


a 


A"' 1 (t-x q) 


(t-x o) 2 


(76) 


Proof: See Appendix lAl 


Given 


M 


at iteration l. according to Lemma [TO] we know 


that |rfc| p (p > 2) is majorized at 


Si) 


[M] by 


Ofe|ffc| 2 + b k \r k \ + a k 


Si) 


- 0 - 1 ) 


Si) 


(77) 


where 


t p - 


a k = 


Si) 




p-i 


(* - r k ] ) 


b k = P 


St) 


p -1 


0- 

2U/j; 


St) 


St) 


(78) 

(79) 


Since the objective decreases at every iteration in the 
MM framework, at the current iteration l, it is sufficient 
to majorize |rfc| p over the set on which the objective is 


smaller, i.e., M" < Efi 


JV-l 


M 


which implies 


Hence we can choose t = 


Vk\ < (Ef=i 1 r% )P y. 

(^2^=1 r k? ) P Then the majorized problem of (l74t 

in this case is given by (ignoring the constant terms) 


N -1 

minimize E {a k \r k \ 2 + b k \r k \) 

k =1 

subject to \x n \ = 1, n = 1,..., TV. 

We can see that the first term Efc=-/ Ofc|nc| 2 in the objective 
is just the weighted ISL metric with weights w k = a k , and 
thus can be further majorized as in Section 1III-B1 Following 
the steps in Section lHI-Bl until ([2(1 (i.e., just the first majoriza¬ 
tion), we know that it is majorized at x(*) by (with constant 
terms ignored) 

x H (r - A max (L)x (i) (x (Z) ) ff ) x, (81) 

where R and L are defined in (ITTI i and (fT4l > with w k = a k . 
For the second term, since it can be shown that b k < 0, we 
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have 


N-l 


N—l 


y> \rk i < Y • 


JO 


fc=l 


fc=l 

N-l 




= Y S Tr(U-fcxx^) 


fe=i 


= Re < x 


if 


2 X 


H 


( N-l 

k=l 

N-l 

E »* 

l)c=1-JV 


„(0 


.(0 


u_ 


„(0 


.(0 


u_ 



(82) 


(83) 


(84) 


(85) 


where = bk,k = 1,..., TV — 1, &o = 0. By adding the 
two majorization functions, i.e.. (|8D and ([85), and defining 


w- k =w k =a k + 


b k 


„(0 




P~ 2 


, k = 1,..., N — 1, 

( 86 ) 


we have the majorized problem of (180) given by 


minimize 


r H 


subject to 


(R-A max (L)x«(x«) ff ) 
= 1, n = 1,... ,N, 


where 


N-l 

E 

k=l-N 


r ( l l- 


(87) 


( 88 ) 


We can see that the problem (187) has the same form as (120) . 
then by following similar steps as in section IIII-BI we can 
perform one more majorization step and get the majorized 
problem 


minimize 

ll x -y|| 2 


X 

(89) 

subject to 

\x n \ 1, Tl 1,..., N, 

y — (Amax 

(L)7V + A u )xW -RxW, 

(90) 


where 


and this time A max (L) should be computed based on weights 
afc in (178) and X u is based on the weights Wk in (186) . As 
in previous cases, we only need to solve ( [89b in closed form 
at every iteration and the overall algorithm is summarized in 
Algorithm[4] Note that, to avoid numerical issue, we have used 
the normalized a* and Wk (i.e., divided by t p ) in Algorithm]]] 
which is equivalent to divide the objective in (180) by t p during 
the derivation. It is also worth noting that the algorithm can 
be accelerated by the scheme described in subsection I V-B I 


VII. Numerical Experiments 

To compare the performance of the proposed MWISL 
algorithm and its variants with existing algorithms and to show 
the potential of proposed algorithms in designing sequences 
for various scenarios, we present some experimental results in 
this section. All experiments were performed on a PC with a 
3.20GHz i5-3470 CPU and 8GB RAM. 


Algorithm 4 Monotonic minimizer for the /: y) -metric of auto¬ 
correlation sidelobes (j> > 2). 


Require: sequence length TV, parameter p > 2 
1: Set l = 0, initialize x(°). 

repeat 


7 

8 
9 

10 

11 

12 : 

13: 

14 

15 


f = F[xW T ,0 lxJV [ 
r=^F"|f | 2 
t = ||r 2 :2v| 


_ l + (p-l)(l^±li) P -p(l^±li) P 1 _ 


° fc _ (Hr,c+llE 

^fc=2 2 ,k = l,...,N -1 


k = 1,. 

-1 

\ L = maxfc {ak(N — k)\k = 1,..., N- 1} 

c = ro [0, wi,..., w N -i,0, u)jv-i, • • •, u>i} T 

(i = Fc 


A„ = \ ( max fi 2 i 


max p 2 i-i) 
l<i<N ' 


_ (l) _ F ,i:xr(A°f) 

y — A 2N{X L N+X a ) 


x ( n +1) = e jarg 

l 4-1 + 1 

until convergence 


n = 1,... ,7V 


,JV-1 


A. Weighted ISL Minimization 

In this subsection, we give an example of applying the 
proposed MWISL and MWISL-Diag algorithms (with and 
without acceleration) to design sequences with low correlation 
sidelobes only at required lags and compare the performance 
with the WeCAN algorithm lfl 8 l . The Matlab code of the 
benchmark algorithm, i.e., WeCAN, was downloaded from the 
websit^H of the book {5l . 

Suppose we want to design a sequence of length N = 
100 and with small correlations only at ri,...,r 2 o and 
7 " 5 i,..., r 70 • To achieve this, we set the weights as 

follows: 

fl, fc € {1,..., 20} U {51,..., 70) (9[) 

[0, otherwise, 

such that only the autocorrelations at the required lags will be 
minimized. Note that in this example there are N — 1 = 99 
degrees of freedom, i.e., the free phases of { x n ln=t (*e 
initial phase does not matter), and our goal is to match 80 
real numbers (i.e., the real and imaginary parts of r \,..., r- 2 o 
and 7*51, • ■ ■, T 70 ). As there are enough degrees of freedom, 
the weighted ISL can be driven to 0 in principle. So in this 
experiment, we will allow enough iterations for the algorithms 
to be run and will not stop until the weighted ISL goes below 
IQ-to r |h e algorithms are initialized by a same randomly 
generated sequence. The evolution curves of the weighted ISL 
with respect to the number of iterations and time are shown in 
Fig- [T|and [2] respectively. Lrom Lig. Q] we can see that all the 
algorithms can drive the weighted ISL to 10 -1 ° when enough 
iterations are allowed, but the proposed algorithms (especially 
the accelerated ones) require far fewer iterations compared 
with the WeCAN algorithm. In addition, we can see that the 
MWISL-Diag algorithms (with and without acceleration) con¬ 
verge a bit faster than the corresponding MWISL algorithms, 

1 http://www.sal.ufi.edu/book/ 
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which means the majorization function proposed in Section 
HVl is somehow better than the one in Section IIII-BI From 
Fig. [2] we can see that in terms of the computational time 
the superiority of the proposed algorithms is more significant, 
more specially the accelerated MWISL and MWISL-Diag 
algorithms take only 0.07 and 0.06 seconds respectively, while 
the WeCAN algorithm takes more than 1000 seconds. It is 
because the proposed MWISL (and MWISL-Diag) algorithms 
require only four FFT operations per iteration, while each 
iteration of WeCAN requires N computations of 2A'-point 
FFTs. The correlation level of the output sequence of the 
accelerated MWISL-Diag algorithm is shown in Fig. [3] where 
the correlation level is defined as 


correlation level = 20 log 10 


n 

ro 


, k = 1 - N, ..., TV - 1. 


We can see in Fig. [3 that the autocorrelation sidelobes are 
suppressed to almost zero (about -160dB) at the required lags. 



Figure 3. Correlation level of the sequence of length N = 100 designed by 
accelerated MWISL-Diag algorithm with weights in ED- 



Figure 1. Evolution of the weighted ISL with respect to the number of 
iterations. 



Figure 2. Evolution of the weighted ISL with respect to time (in seconds). 
The plot within the time interval [0,10] second is zoomed in and shown in 
the upper right comer. 


11. PSL Minimization 

In this subsection, we test the performance of the proposed 
Algorithm [3 in Section [VI] in minimizing the peak sidelobe 
level (PSL) of the autocorrelation sidelobes, which is of 
particular interest. To apply the algorihtm, we need to choose 
the parameter p. To examine the effect of the parameter p, we 
first apply the accelerated version of Algorithm |3 (denoted as 
MM-PSL) with four different p values, i.e., p = 10,100,1000 
and 10000, to design a sequence of length N = 400. Frank 
sequences ED are used to initialize the algorithm, which 
are known to be sequences with good autocorrelation. More 
specifically, Frank sequences are defined for lengths that are 
perfect squares and the Frank sequence of length N = M 2 is 
given by 

XnM+k+i = e j27rnfe / M , n, k = 0,1,..., M — 1. (92) 

For all p values, we stop the algorithm after 5 x 10 4 iterations 
and the evolution curves of the PSL are shown in Fig. ©. 
From the figure, we can see that smaller p values lead to faster 
convergence. However, if p is too small, it may not decrease 
the PSL at a later stage, as we can see that p = 100 finally 
gives smaller PSL compared with p = 10. It may be explained 
by the fact that £ p -norm with larger p values approximates 
the Loc-norm better. So in practice, gradually increasing the p 
value is probably a better approach. 

In the second experiment, we consider both an increasing 
scheme of p (denoted as MM-PSL-adaptive) and the fixed p 
scheme with p = 100. For the increasing scheme, we apply the 
MM-PSL algorithm with increasing p values 2, 2 2 ,..., 2 13 . 
For each p value, the stopping criterion was chosen to be 
|obj(x( fc+1 )) — obj(x( fe ))| /obj(x( fc )) < 10 -5 /p, with obj(x) 
being the objective in (T73l >. and the maximum allowed number 
of iterations was set to be 5 x 10 3 . For p = 2, the algorithm 
is initialized by the Frank sequence and for larger p values, 
it is initialized by the solution obtained at the previous p. 
For the fixed p scheme, the stopping criterion was chosen 
to be |obj(x( fe+1 )) — obj(x( fe ))| /obj(x^) < 10 -1 °, and the 
maximum allowed number of iterations was 2 x 10 5 . In this 
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Figure 4. The evolution curves of the peak sidelobe level (PSL). 


case, in addition to the Frank sequence, the Golomb sequence 
ESI was also used as the initial sequence, which is also known 
for its good autocorrelation properties. In contrast to Frank 
sequences, Golomb sequences are defined for any positive 
integer and a Golomb sequence {x n }^ =1 of length N is given 
by 

x n = e i 7r ( n - 1 W JV ) n = 1,... ,N. (93) 

The two schemes are applied to design sequences of the 
following lengths: N = 5 2 , 7 2 ,10 2 , 20 2 ,30 2 , 50 2 , 70 2 ,100 2 , 
and the PSL’s of the resulting sequences are shown in Fig. 
[3 From the figure, we can see that for all lengths, the MM- 
PSL(G) and MM-PSL(F) sequences give nearly the same PSL; 
both are much smaller than the PSL of Golomb and Frank se¬ 
quences, while a bit larger than the PSL of MM-PSL-adaptive 
sequences. For example, when N = 10 4 , the PSL values of the 
MM-PSL(F) and MM-PSL-adaptive sequences are 4.36 and 
3.48, while the PSL values of Golomb and Frank sequences 
are 48.03 and 31.84, respectively. The correlation level of 
the Golomb, Frank and the MM-PSL-adaptive sequences are 
shown in Fig. [6] We can notice that the autocorrelation 
sidelobes of the Golomb and Frank sequences are relatively 
large for k close to 0 and N — 1, while the MM-PSL-adaptive 
sequence has much more uniform autocorrelation sidelobes 
across all lags. 


VIII. Conclusion 

We have developed two efficient algorithms for the min¬ 
imization of the weighted integrated sidelobe level (WISL) 
metric of unit-modulus sequences. The proposed algorithms 
are derived based on applying two successive majorization 
steps and we have proved that they will converge to a sta¬ 
tionary point of the original WISL minimization problem. By 
performing one more majorization step in the derivations, we 
have extended the proposed algorithms to tackle the problem 
of minimizing the f p -norm of the autocorrelation sidelobes. 
All the algorithms can be implemented by means of FFT 
operations and thus are computationally very efficient in 


Figure 5. Peak sidelobe level (PSL) versus sequence length. MM-PSL(G) 
and MM-PSL(F) denote the MM-PSL algorithm initialized by Golomb and 
Frank seauences. resnectivelv. 
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Figure 6. Correlation level of the Golomb, Frank and MM-PSL-adaptive 
sequences of length N = 10 4 . 


practice. An acceleration scheme that can be used to further 
speed up the proposed algorithms has also been considered. By 
some numerical examples, we have shown that the proposed 
WISL minimization algorithms can generate sequences with 
virtually zero autocorrelation sidelobes in some specified lag 
intervals with much lower computational cost compared with 
the state-of-the-art. It has also been observed that the proposed 
/' r ,-metric minimization algorithm can produce long sequences 
with much more uniform autocorrelation sidelobes and much 
smaller PSL compared with Frank and Golomb sequences, 
which are known for their good autocorrelation properties. 
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Appendix A 
Proof of Lemma[Tq1 

Proof: For any given xq € [0, t), let us consider a 
quadratic function of the following form 

9{x\x 0 ) = f{x o) + /'(x 0 )(x - x Q ) + a(x - x 0 ) 2 , (94) 


where a > 0. It is easy to check that f(x o) = g{xo\xo). So to 
make g(x|xo) be a majorization function of f(x) at xq over 
the interval [ 0 , t], we need to further have f{x) < g(x |xo) for 
all x £ [0, t\, x ^ xq. Equivalently, we must have 

^ f(x) - f(x 0 ) - /'(x 0 )(x - x 0 ) 

“ (x- xn ) 2 { 

yas o/(j j 

for all x € [0, t], x ^ xq. Let us define the function 

^ f{x) - f{x 0 )~ f'{x o)(x-x 0 ) 

A(x\x 0 ) = - 7 -72- ( 96 ) 

(x - x 0 p 

for all x ^ Xo- The derivative of A(x|xo) is given by 

A >r x \ X0 \ = PW + ~ 2 (/( x ) ~ f( x °))/( x ~ x o) 

0 (x x 0 ) 2 

Since f(x) = px p ~ 1 is convex on [0, t] when p > 2, we have 

= /'/'(*„ + T(x - xa)) ir 

X-Xq J o 

< f (/'(x 0 )+ t(/'(x) -/'(x 0 )))dT 
Jo 

= 7; if (x) + f (x o)), 

which implies A'{x\xq) > 0 for all x G [0,f],x 7 ^ xo- Thus, 
A(x|xo) is increasing on the interval [ 0 , t\ and the maximum 
is achieved at x = f. Then the smallest a we may choose is 


a = max A(x|xo) 

x£[0,t],x^xo 

_ t p - Xg -pxg -1 (f - Xq) ( ' 97 ' ) 

(f - Xo ) 2 

By substituting a into g(x\xo) in (l94t and appropriately 
rearranging terms, we can obtain the function in ( T75l i. ■ 
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